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Abstract 

A method of fundamental solutions has been used to study adiabatic transition ampli- 
tudes in two energy level systems for a class of Hamiltonians allowing some simplifications 
of Stokes graphs corresponding to such transitions. It has been shown that for simplest 
such cases the amplitudes take the Nikitin - Umanskii form but for more complicated 
ones they are formed by a sum of terms strictly related to a structure of Stokes graph 
corresponding to such cases. This paper corrects our previous one [Phys. Rev. A, 63 
052101 (2001)] and its results are in a full agreement with the ones of Joyc, Mileti and 
Pfister [Phys. Rev. A, 44 4280 (1991)]. 

PACS number(s): 03.65.-W , 03.65. Sq , 02.30.Lt , 02.30.Mv 

Key Words: fundamental solutions, semiclassical expansion, JWKB approximations, exponential 
asymptotics, adiabatic limit 

1 Introduction 

In our previous paper [J we have applied a formalism of fundamental solutions to obtain 
formulae for adiabatic transition amplitudes in two level energy systems. The corresponding 
formalism has been developed under quite general assumptions about a nature of Hamil- 
tonians perturbing a system adiabatically. Unfortunately, in its applications to particular 
examples considered in the section V and the furthers of the paper mentioned we have made 
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a fatal error in detailed calculations of the corresponding transition amplitudes. As a conse- 
quence of this we have also drawn in this paper erroneous conclusions which followed from the 
obtained erroneous formulae. In the present paper we would like to correct the corresponding 
calculations as well as to draw correct conclusions. 

However in order to avoid a permanent referring to the material presented and discussed 
in the first four sections of the paper [I] and to make our present paper selfsufficient and 
selfconsistent we shall repeat below to large extent the contents of these sections. Therefore 
we shall start with reminding shortly main reasons for studying transitions in two energy 
level systems. 

First of all such systems provide us with the simplest models to investigate transition 
amplitudes between different energy levels by different approaches [2]. On the other side 
these systems play an important role in experimental investigations of basic principles of 
quantum mechanics j^j. Recently a lot of effort has been devoted to obtain more rigorous 
results on the adiabatic limit of transition amplitudes for these systems jU E3 HI E] ■ In 
particular in a series of recent papers Joye et al @J |SJ H3 HI have studied this problem 
by the Hilbert space methods. Such two energy level systems are formally equivalent to a 
one-half spin system put into time dependent magnetic field. However good approximate 
results and the more so the exact ones are difficult to obtain for such systems even for simple 
time evolutions of the effective 'magnetic' field. Therefore each opportunity of improving this 
situation is worth trying. A treatment of the problem by a method of fundamental solutions 
(so fruitful in its application to stationary problems of 1-dim Schrodinger equation ) 
is of first importance, the more so that to our knowledge, the method was not used so far 
to this goal. A possibility of application of the method is related to the fact that a linear 
system of first order differential equations describing time evolution of transition amplitudes 
can always be transformed into a system of decoupled second order equations having a form 
of the stationary Schrodinger equation, one for each amplitude. This allows us to apply all 
advantages of the fundamental solution method 11 j. The only obstacle related with this 
approach is a complexity of effective 'potentials' which appear in the final system of the 
Schrodinger-type equations. 

The paper is organized as follows. 

In the next section the problem of transitions in two energy level systems is stated and 
corresponding assumptions about the effective 'magnetic field' are formulated. A linear sys- 
tem of two differential equations for the transition amplitudes is rewritten in a form of two 
decoupled equations of the Schrodinger type. 

In Sec. 3 properties of the fundamental solution method are recalled. 

In Sec. 4 some subtleties of the application of the fundamental solution method to the 
problems considered in the paper are discussed. 

In Sec. 5 a class of Hamiltonians with so called NED property is distinguished for further 
considerations. 

In Sec. 6 an exact form of a transition amplitude for the NED systems is obtained and its 
adiabatic limit is found. 

In Sec. 7 two examples of the NED systems are considered and the Nikitin - Umanskii 
formula is reconstructed. 

Finally in Sec. 8 we discuss our results stressing their coincidence with the corresponding 
ones of Joye et al 0. 
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2 Adiabatic transitions in two energy level systems 



First let us remind that, in general, any two energy level system is formally equivalent 
to a one-half spin system put into an external magnetic field B(i). Its Hamiltonian H{t) is 
given then by H{t) = |//B(i) • a , where a = (a x ,a y ,a z ) are Pauli's matrices so that two 
energy levels E±(t) of H{t) are given by E±(t) = ±%B(t) where B(t) = V /B 2 (t). 

When adiabatic transitions between the two energy levels E±(t) are considered then the 
following properties of the field B(t) are typically assumed: 

1° B(t) is real being defined for the real t, — oo < t < +00; 

2° B(i) can be continued analytically off the real values of t as a meromorphic function 
defined on some t-Riemann surface Re- A sheet of Rg from which B(i) is originally continued 
is called physical; 

3° On the physical sheet B(t) is analytic in an infinite strip £ = {t : |3f£| < <5, <5 > 0}, 
without roots in the strip and achieves there finite limits for 3f?t = ±00 , i.e. B(!Ri = ±00) = 
B± 7^ in the strip; 

4° -B(i) = ^/B 2 (i) is a ramified function of t on Rg with square root branch points 
coinciding with crossing points of the two energy levels E±(t); 

The field B(t) depends additionally on a parameter T(> 0) i.e. B(i) = B(t, T) which 
introduces a "natural" scale of time to the system, so that its time evolution is expressed 
most naturally in units of T. If T is small in comparison with the actual period of the process 
considered then the latter is "fast" or "sudden". If, however, T is large in this comparison 
then the process is " slow" or " adiabatic" . 

In the adiabatic process of the system the following is assumed about B(t,T): 

5° A dependence of B(t,T) on T is such that a rescaled field B(sT, T) has the following 
asymptotic behavior for T — > +00 

B(sT, T) ~ B (s) + |bx( S ) + ^B 2 (s) + ■■■ (1) 

while its s-Riemann surface Rb(T') approaches 'smoothly' the topological structure of the 
Riemann surface corresponding to the first term Bo(s) of the expansion (^Q). 

6° With respect to its dependence on s the field B (s) satisfies properties 1° — 4° above 
with substitutions t — ► s and B(s) — > Bq(s). 

7° For purposes of this paper we shall assume also an algebraic dependence of B(sT, T) 
on s so that its asymptotic behaviour in the strip S as s — ► ±00 is the following: 

B( S T,T)~B ± + 5| + ^ + . + + 

(2) 

1 

- < a\ < a>2 < • • • < afe < • • • 

where «i, . . . , a^j are assumed to be rational. 

The next assumption which validity becomes clear in Sec. 6 needs to formulate a notion 
of Stokes lines for the function Bq(s). These are lines which starts at roots of Bq(s) and are 
governed by conditions (i Bo(a)da^j = where Sk, k = 1,..., are roots of Bq(s). We 
shall assume the following about roots of Bq(s), its Stokes lines and about components of the 
limiting field Bq(s). 

8° 
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bf a) Roots of B 2 (sT,T) and B$(s) are simple. 

b) A boundary of the central strip £ = {t : |9f£| < 5, 5 > 0} in which B(sT,T) is 
holomorphic and with no roots inside it consists of two infinite Stokes lines which become 
complex conjugated with each other in the adiabatic limit. 

c) Each component of B(sT, T) and Bq(s) is holomorphic and nonvanishing at roots of 
B 2 (sT,T) and B 2 (s) respectively. 

The time-dependent Schrodinger equation corresponding to H{t) takes a form 



i d^{s,T) 
T ds 



~»B(sT,T) ■ <tV(s,T) 



(3) 



The adiabatic regime of evolution of the wave function ^f(s, T) corresponds now to taking 
a limit T — > +oo in (j3J). 

The main problem of the adiabatic limit in the considered case is to find in this limit 
the transition amplitude between the two energy levels of the system for s — ► +oo under the 
assumptions that \]/(— oo,T) coincides with one of the two possible eigenstates 5 r ±(— oo, T) 
of H(—oo) (= H(+oo)) (corresponding to E±(— oo) (= E±(+oo)) and that there is no level 
crossing for real t i.e. liminf B(t) > e > 0. Known approximate solutions of this problem 

— oo<t<+oo 

are that of Landau ^5] and Zener ^3] in a form of the so called Landau-Zener formula and 
that of Dykhne ^3] who have shown that such an amplitude should be exponentially small 
in the limit T — ► +oo. In the next sections we shall show how to get an exact (i.e. not 
approximate) result for this amplitude as well as its adiabatic limit with the help of the 
fundamental solutions . 

A typical way of proceeding when the adiabatic limit is investigated is using eigenvectors 
^±(s,T) of H(sT,T) satisfying (vJ/ij^i) = 0. Then, such eigenvectors ^±(s,T) can be 
chosen as the following ones 



COS y 



sinfe* 9 



¥_(s,T) = e 



-i Jq c/>cos 2 ^da 



Sm 2 



cos f e i? 



(4) 



where and <f> are respective polar and azimuthal angles of the vector B(t, T) and dots over 
different quantities mean derivatives with respect to s-variable. 
The wave function Vl/(s,T) can now be represented as 



9(a, T) = o+(a, T)e" iT 'So E+ ^' T)d ^ + (s, T) + a_(s, T)e~ iT So E -^^_( a , T) 



(5) 



The Schrodinger equation (j3J) can be rewritten in terms of the coefficients a±(s, T) as the 
following linear system of two equations 



a + (s,T) = -c*(s,T)e i So UJ ^ T ^ a .(s 1 T) 
d_(a,T) = c(s, r)e _i /o " ( «' T) ^a+(s, T) 



(6) 



where 



Bx(BxB) 



B x B 



c( s ,T) = - + ^ s in9 = --^ ' J Z + 

2 2 2 B\MTBl *BjB*^m 
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(7) 



uj(s,T) =T(E+- E-) - (pcosQ = fiTB 



B x B 



B Bl + Bl 



According to our assumptions we are looking for a solution to the system (JJJJ) satisfying 
the following initial conditions a+(— oo,T) = 1 and a_(— oo,T) = and under this condition 
we are interested in the limits lim^+oo a_(s,T) and limT^ +00 a_(+oo,T). 

It is seen from (JJJ) however that a dependence of the coefficients c and wons and T can 
be quite complicated since the corresponding dependence of these coefficients on the B-field 
components is sufficiently complicated. We can simplify the latter dependence by a suitable 
unitary transformation applied to ^(s, T) in the formula © defined by the following operator 



U = e §VT f° B z (sT,T)dsa z 



ii^T J°B z (sT,T)ds 







\iliT J°B z {sT,T)ds 



(8) 



For the new amplitudes 



ai(s,T) 
a 2 (s,T) 



we get 



/ rp\ -i P (0+/xT_B) sin 2 %da © t m\ , 

ai(s,T) = e J o yv p ' 2 cos 0+ ^ S) j 1 ) + 

-l rU-flTB) COS 2 f dcr .© / T ; 

e Jo vr r ; 2 sm — a-(s,l) 



(9) 



a 2 ( s , T) = e < cos2 1 S in |a+(,, T) - 

e i j>+^)sin 2 f *,+i*(0) CQS ©^ T) 

The transformation © does not change the form of Eq.(jHJ) changing only the correspond- 
ing functions c and u. Namely we have 



di (s,T) = -Tcl(s,T)e iT So ^ T ^a 2 (s,T) 



a 2 (s,T)=Tc 1 (s,T)e~^o^ T)d ^a 1 (s,T) (10) 

where 

Cl = - ~fj,B sin 6e i<? \ wi = /iBcosG (11) 

It is worth to note that the form © of the considered transformation provides us imme- 
diately with the asymptotic forms of the amplitudes ai(s,T) and a 2 (s,T) for s — ► —00 since 
the amplitudes a±(s, T) can take arbitrary values a±{— 00, T) in this limit satisfying only the 
condition |a+(— 00, T)| 2 + |a_(— oo,T)| 2 = 1. Namely we have simply in this limit 

/ rTi\ -i ["(<i>+uTB) sin 2 %da ® / rri\ , 

ai(s,T)~e J o ^ p ; 2 cos — -a+(— oo,T) + 
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-i f°(<i>-uTB) cos 2 % da ■ © / rp\ 

e Jo ^ ^ ; 2 sin — a_(— oo, i ) 



2_( — OO, 'J'j 

(12) 



I m\ i f S '(i-uTB) cos 2 #dcr+j<A(0) • © / rp\ 

a 2 (s,lj ~ e Jo ^ p ; 2 T sm a+(— 00, T) — 

t ru+uTB) sin 2 S(ij+i<A(0) © / T \ 
e J y 2 COS — (X_ ( — OO, 1 ) 

However since we are going to consider the case a + (— 00, T) = 1 and a_ (— 00, T) = then 
for this case we get for s — > — 00 

/ -i f 3 U+uTB) sin 2 %dcr © 

ai(s,j)~e Jo ; 2 cos — 

I rr\ i f 3 (<j)- uTB) cos 2 ^-da+icbtO) • © /-, \ 

a 2 (s, i ) ~ e J o vv p ; 2 T sin — (13) 

We can express further the amplitude a_(s,T), in which we are interested, by the a± 2 
ones inverting the transformation © to get 

e 



I rr\ i f„ S ii-nTB) cos 2 #d<r . " . . 

.(s,T) = e J o^ p ' 2 sm —a\{s, T) + 
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-i f s (<j>+uTB) sin 2 §dcr-i<i(0) © / A \ 

e Jo ' 2 cos _ a 2 {s,T) (14) 

Moreover we can always assume that a limiting value of the field B(sT, T) for s — ► +00 
coincides with its z-component to be B + = B(+oo,T) = (0,0,73+), so that 6(+oo,T) = 0. 
Therefore in the limit considered we get from (|14l) 

a„(+oo,T) = lim e^Io^^ sin 2 1*^(0)^ T) (15) 

s— >+oo 

Consequently it is the amplitude 02(5, T) for which the above limit we have to consider. 
The system (jlOj) can be rewritten further as the following linear system of second order 
equations 



c* \ 

-r + iTwi ) di + IciPai = 
c\ J 



ai 

a 2 - - iTuj^j d 2 + |ci| 2 a 2 = (16) 

where the amplitudes ai j2 decouple from each other being however still related by ©• 
By the following transformations 



a l {s,T)=e ZJ0 \^ ) h(s,T) 

a 2 (s,T) = J I °&- iT " 1 )%(s,T) (17) 
we bring the equations to Schrodinger types 

6 1)2 (s, T) + TV, 2 (s, T)b h2 (s, T) = (18) 
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where 



Q2(s,T) 



1 



ci 



iTtJi +|ci| + 



1 



4T 2 Vci 

while (for real s and T) we have 

gi(s,T) = ql{s,T) 
A dependence of the function q2(s,T) on T is given by 



ci 



2T 2 Vci 



92(s,r) = -n B - I — I + 



1 



2T B« 



ci 



2T 2 



1 fty 



2 Vci 



(19) 



(20) 



(21) 



where a dependence of B,B z ,c\ on T in (|21|) is also anticipated. By 1)20(1 we get a corre- 
sponding dependence of qx(s,T) on T. 

The equations 1)18(1 are now basic for our further analysis since their form is just of the 
stationary 1-D Schrodinger equation. 

Taking into account (JJ) and Q it is easy to check that the last formula provides us with 
the following type of asymptotic behavior of 91.2 (s> T) for large T: 



q l>2 (s,T) = + + ^$00 + • • • 



(22) 



Therefore the above form of dependence of qi^(s,T) on T permits us to apply to the 
considered case the method of fundamental solutions. For this reason we shall start the next 
section with a review of basic principles of the method suitably adapted to the considered 
case. 



3 Fundamental solutions and their properties 

Consider first qi t 2(s, T) as functions of s. They are defined completely by an s-dependence 
of field B(Ts, T). According to our assumptions, the latter is meromorphic on some Riemann 
surface Rb(T). However, by ()2*T)) . qi^(s,T) are algebraic functions of B, B and B and, there- 
fore, they are also meromorphic functions of s defined again on some other Riemann surfaces 
Hi 2 determined by these algebraic dependencies. As it follows from ()21j) topological struc- 
tures of Ri,2 can be quite complicated. However, in what follows, we are interested in the 
adiabatic limit T — > +00 by which structures of Ri 2 should be determined for T large enough 
basically by the first term q^\{s) of the expansion (|22j). In consequence, by ()2*T)) . it should 
be determined by /xB^(s) i.e. by the first term of the expansion (^Q). Structures of Ri,2 can 
turn out to be much simpler in this limit. 

Despite this supposed complexity of 13^2 (s,T) and of their Riemann surfaces we shall 
introduce and discuss the fundamental solutions to the equations 1)18(1 without simplifications. 
We shall do it for the q2(s,T) case). An extension of the discussion to the q\(s,T) case will 
be obvious. 

A standard way of introducing the fundamental solutions is a construction of a Stokes 
graph EH EU related to a given (72 (s, T). Such a construction, according to Frdman and 
Froman |5] and Fedoriuk 10 , can be performed in the following way jllj . 
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Let Z denote a set of all the points of R2 at which q2(s, T) has its single or double poles. 
Let 5(x) be a meromorphic function on R2, the unique singularities of which are double poles 
at the points collected by Z with coefficients at all the poles equal to 1/4 each. (In a case 
when R2 is simply a complex plain the latter function can be constructed in general with the 
help of the Mittag-Lefner theorem ^5] . But for a case of branched R2 the general procedure 
is unknown to us). Consider now a function 

q 2 (s,T)=q 2 (s,T) + ^5(s) (23) 

The presence and the role of the <5-term in 1)23|) are explained below. This term contributes 
to (|23|) if and only when the corresponding 'potential' function q2(s,T) contains simple or 
second order poles. (Otherwise the corresponding 5-term is put to zero). It is called the 
Langer term ^2 IS] • 

The Stokes graph corresponding to the function q%(s,T) consists now of Stokes lines 
emerging from roots (turning points) of q2(s,T). Stokes lines satisfy one of the following 
equations: 

9 f Jm,T)dH = (24) 

J Si 

with Si being a root of (72(5, T). We shall assume further a generic situation when all roots 
Si are simple. 

Stokes lines which are not closed end at these points of R2 (i.e. have the latter points 
as their boundaries) for which the action integral in Q24|) becomes infinite. Of course such 
points are singular for q~2(s, T) and they can be its finite poles or its poles lying at an infinity. 

Each such a singularity of q~2(s,T) defines a domain called a sector. This is the con- 
nected domain of R2 bounded by Stokes lines and Zk itself. The latter is also a boundary for 
the Stokes lines or being an isolated boundary point of the sector (as it is in the case of the 
second order pole). 

In each sector the LHS in (|24jl is only positive or only negative. 

Consider now equation (|18|) for 62(5, T). Following Froman and Froman in each sector Sk 
having a singular point Zk at its boundary one can define a solution of the form: 

M*> T) = f 2 l {s, T)-e° iTW ^ X2 ,k{s, T) k = 1, 2, . . . (25) 

where 

X 2,k{s, r) = 1 + E (-^rT f d6 f 1 d ^ 2 ■ ■ ■ f" 1 ^» n (^)n(6) • • • * 

n>l ' z h z h z k 

(26) 

'l _ e -2aiT{W{s)-W{^))^ ^ _ e -arirCW(£i)-W(&))^ . . . (l _ e -2viT(W(t n -i)-W(S n ))^ 



with 



and 



^ S(s) lqo(s,T) 5 q' 2 2 (s,T) 
qUs,T) %~|( S ,T) lb qUs,T) 



W(s,T)= Jq(£,T)dZ (28) 



8 



where s, is a root of q(s,T) lying at the boundary of S^. 

In (|25|) and (|27j) a sign of cr (=±1) and an integration path are chosen in such a way to 
have: 



a<3(W(Z 3 )-W(Z j+1 ))<0 



(29) 



for any ordered pair of integration variables (with £o = s). Such an integration path is then 
called canonical. Of course, the condition (|29|) means that 62^(5, T) vanishes in its sector 
when s — > z^ along the canonical path. The Langer <5-term appearing in (|23|) and (|27|) is 
necessary to ensure all the integrals in ()27|) to converge when is a first or a second order 
pole of <j2(s,T) or when the solutions (|23)) are to be continued to such poles. As it follows 
from (|27|) each such pole Zk demands a contribution to 5(s) of the form (2(s — what 
has been already assumed in the corresponding construction of S(s). 

It is now necessary to mention the main property of the fundamental solution method 
which is that analytic continuations of fundamental solutions along canonical paths ensures 
an immediate pass to adiabatic limit on every stage of calculations performed with their use. 
This property can be always utilized if all zeros of ^-functions are simple and distributions of 
their zeros and poles are discrete i.e. there are no accomodation points for these singularities. 
We shall assume in the remainder that the two level energy systems we are going to consider 
will satisfy the last conditions. 

4 The adiabatic limit in the fundamental solution approach 

Consider now the consequences of taking the large-T limit for the above description. 
We assume that for a given ij2(s,T) and its Riemann surface R2 the corresponding Stokes 
graph G2 is drawn. It is drawn, of course, on the Riemann surface corresponding to 



First let us notice that singular points of q2(s,T) such as its branch points and poles 
depend in general on T. For both kinds of these singularities this also means a dependence 
on T of jumps of (72 (s> T) on its cuts as well as a T-dependence of coefficients of its poles. 

According to the property 5° of the magnetic field B (see Sec. 2) we can expect that a 
singular structure of <j2(s,T), i.e. positions of its roots and poles, as well as cut jumps and 
pole coefficients, change smoothly in this limit to their final positions and values respectively. 
This limit structure is defined by a singularity structure of q% (s,T) (see expansion (I22J1 ). 
Therefore, both a topology of VR2 and the associated Stokes graph G2 changes accordingly 



to coincide eventually with the Riemann surface y R 2 and with the Stokes graph G 2 



corresponding to y q\ (s, T) . This limiting structure can be achieved in the following ways: 



a. some of branch points and poles of q2(s, T) escape to infinities of R2; 

b. some of branch points and poles of q+(s,T) approach the respective singularities of 

Q2°\s,n 

c. some of branch points and poles of (j2(s,T) disappear because their respective jumps 
and coefficients vanish in the limit T — > +00. 

Being more specific we expect that for T large enough a set S2 of all singular points 
of <j2(s,T) (i.e. containing all its branch points and poles) consists of three well separated 
subsets S 2 ™^, S2 an and S^*™- The points of S 2 n ^ run to infinities of R2 when T — > +00. Those 
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of S v 2 an 



disappear in this limit while those of S 2 
singular points of q2(s,T) . 



fin 



coincide in this limit with the set sl°^ of 



Let us remove the points contained in S 2 



m f U S 2 °™ from the Riemann surface R2 , i.e. 



let us consider these points as regular for q2(s,T). Then R2 will transform into R 2 m - a 
Riemann surface which singular points coincide with those of the set S 2 

Together with the previous operation let us remove from \fRv also the Stokes lines gener- 
ated by the points of S 2 n/ U S% an so that the remaining Stokes lines can be uniquely continued 



fin 



to form the Stokes graph G 2 m generated by the set S 2 m . It is clear that the graph G 
coincides with G 2 °^ in the limit T — > +00. 

The above two operations will be called the adiabatic limit reduction or simply the re- 
duction operation. 

As we have mentioned earlier there is a set of sectors and a corresponding set of funda- 
mental solutions defined in them associated with the graph G2. By the reduction operation 
both sets can be reduced i.e. under this operation some sectors of G2 transform into cor- 
responding sectors of G 2 whereas the others disappear. Obviously, the latter sectors are 
those which disappear when the limit T — > +00 is taken. 

A good illustration for the above discussion can be an example considered in Sec. 7, namely 
the Nikitin model of the atom-atom scattering, for which the corresponding rescaled B-field is 

the following B(sT,T) = ( (b 2 



fin 



,fin 



0, 1 I We can write for this case the corresponding 



Schrodinger equation (|18j) using the amplitudes a±(s, T) for which the respective q±(s, T) are 
following 



2T 2 



q±(s,T)= — 




2s 2 + b 2 (b 2 + s 2 ) 
s 2 (b 2 + s 2 ) 



6s(b 2 + s 



2\2 



1 



2T I 1 + (b 2 + s 2 ) 3 b 2 + s 2 



(30) 



2 T (l + (6 2 + S 2 ) 3 )2 (&2 + s 2)f 

3 4(6 2 + ^ 2 ) 4 (g 2 - b 2 ) - 4(b 2 + s 2 ){b 2 + 5s 2 ) + 3s 2 (b 2 + s 2 ) 

2 (l + ( & 2 + s 2)3)2 



Equations (f30|> show that in the limit T — > +00 the Stokes graph for the considered 
problem is determined by the function 



(Ae) 



1 + 



1 



{b 2 + s 2 )\ 



(31) 



The graph is shown on Fig.l. 

While each of q±(s,T) has 40 roots, five branch points at s = ±ib and at s = Sk = 

(2k + l)iri \ 2 

3 — b \ , k = 1,2,3, as well as two poles at s = 0, there are only six roots at 

s = Sk, k = 1, 2, 3 and only two poles at s = ±ib for q(°\s, T). 

The functions q±(s,T) are determined on two sheeted Riemann surfaces R-t respectively 
with the branch points at s = ±ib and at s = Sk, k = 1, 2, 3 and with 40 roots distributed into 
halves on each sheet of the surfaces. Therefore the Riemann surfaces \/R± corresponding 
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to \/q±(s, T) are four-sheeted with these 40 roots being square root branch points on them. 
When T — * +00 only six of these branch points survive coinciding with the six roots of 
q(°\s,T) at s = =tsfc, k=l,2,3 whereas R-t transforms into the complex s-plane since the 
branch points of q±(s,T) at s = ±ib disappear, being transformed into the second order 
poles of g (0) (s,^). It is easy to check however that for finite but large T these six roots of 
g(°)(s, T) are each split initially into two as. The split is the result of the square root branch 
points at s = ±ib to which the recovering of the finite T transforms the poles of q(°\s, T) at 
the same points. The two copies of each of these six roots lie of course on different sheets of 
H±. Next, each of these 12 roots is still split into three by the same reason of finiteness of 
T. In this way, on each of the two sheets of R-t there are 36 roots grouped by three around 
their limit s = is^, k = 1, 2, 3 achieved for T — > +00. 

The remaining four roots of q±(s,T) are displaced in two pairs, one pair on each sheet of 
R-t, close to the points s = at which the second order poles of q±(s, T) are localized. When 
T — * +00 the roots in each pair collapse into s = multiplying the corresponding second 
order poles and thus causing mutual cancellations of the latter and themselves in this limit. 

Now let us focus our attention on the Stokes graph G_ generated by g_(s,T) on the 
first sheet of R_ as well as on the remaining ones. It looks as in Fig. 2. (The Stokes graph 
G + corresponding to q + (s,T) can be obtained from G_ by complex conjugation of the 
latter.) On the figure the wavy lines denote the cuts corresponding to the branch points 
of the fundamental solutions defined on R_. The sheet on Fig. 2 cut along the wavy lines 
defines a domain where all the fundamental solutions 6_ 5 i(s,T), b_ 2(s, T) defined in the 
corresponding sectors Si, S% (shown in the figure) are holomorphic. 

According to our earlier description of the behavior of the Riemann surface y/R+ when 
T — > +00 the set S!?^ corresponding to the considered case is empty, S!i an contains four 
points at s = on each of the four sheets of ^R_ (these four points correspond to the second 
order poles of q~(s, T)) and the four branch points close to s = 0, while S^ n contains all the 
remaining singular points of y/q_(s, T). 



5 Systems which are not essentially different from their adia- 
batic limits (NED systems) 

The last example considered above shows us also that by changing the amplitude repre- 
sentation to the ai t 2 ones we get much simpler s-dependence for the corresponding functions 
312 defining Eas. (|18|) and for their adiabatic limit q^\- Namely for the interesting us ampli- 
tude 02 we have 

( r) _ h A \2 (-> 1 A 3* Ag f 3 2b 2 + s 2 

Q2{s, ) 4 ^ej + (62 + s2)3 j 2T b 2 + s 2 AT 2 (b 2 + s 2 ) 2 

(32) 

It is seen from (|32[) that both the functions q2(s,T) and q^is) have the same Riemann 
surfaces, namely the simple complex plain on which they have poles in exactly the same 
points. They differ only by positions of their zeros the latter being in a mutual one-to-one 
correspondence so that each zero of q^ (s) is an adiabatic limit of the corresponding zero of 
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q2(s, T). Therefore the Stokes graphs corresponding to both these functions are topologically 
equivalent having the forms of Fig.l 

As a consequence of this an application of the fundamental solution method to the cases of 
Eqs.(|18|) with respective q2(s,T) and #2 ( s ) functions gives the same adiabatic limit for both 
these cases. We shall describe such a situation as corresponding to a system which do not 
differ essentially from its adiabatic limit and we shall call such a system the not-essentially- 
different one (the NED-system). 

It follows from the above discussion that the NED property is not an immanent one for 
a system but can be achieved by choosing a suitable amplitude representation for a system. 




6 Transition amplitudes for NED systems 

As it follows from the discussion of the previous section systems with the NED properties 
allow us for as easy canonical continuations of fundamental solutions of interests as they are 
for their corresponding adiabatically reduced forms. Therefore for such systems we can 
consider them applying an exact procedure or using the adiabatical limit for the latter to get 
correct results for adiabatical limit transition amplitudes. 

We shall apply the procedure of canonical continuation of fundamental solutions to the 
amplitude a 2 (s,T). First we have to express this amplitude by the fundamental solutions 
and to satisfy the second of the conditions (|T3|). Canonically continued to —00 with simple 
results of such continuations are the solutions &2,i( s > T)and b2^(s,T) corresponding to the 
sectors S± and Sj respectively shown in Fig.3 representing a Stokes graph corresponding to 
a general NED system. We have 

a 2 (s,T)=e J ° 2 ^ (A6 2>I ( 8 ,r) + D6 2)1 ( fl> T)) (33) 

where 

M*,T) = «2hs,T)e T ^ V ^ F)d ' X 2As,T) (34) 
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and where we have assumed the positive real value of yq^s, T) on the physical sheet. The 
lower integration limits in (|34j) are respective zeros of q2(s,T) shown in Fig. 3. 




Fig. 3 The Stokes graph corresponding to q2(s,T) of a NED system 
Taking into account that 

ci B 



a B 

B ? B 



+ 9 cot G + it 



■ cos © — 6 sin 
B B 

we get in the limits s — ► ±oo along the real axis 

iTJq 2 (s,T) ~ ~ifj,TB+~^- ( 5* - £l ] = Ji/iTfi - ^— ^- - J#cos9 
v 2 2 B \B Z c\ 2 2 sin B 2 



so that 

2 Ui 
1 (&L 

in both the limits. 



■ . . 2 e is e e 

iu\ ] + iT^fcfr ~ («0 + ^ / uT J B) sin — + - — - — tan — 



, . 2 e IB e e 

ia>i ) - ~ (up - ifiTB) cos — + - — + — cot — 



(35) 



(36) 



(37) 
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It is now easy to show that the Eas. (|33|) -(|37 |) provide us with the following asymptotic 
form of the amplitude a2(s, T) for s — ► — oo 



a 2 (s,T) 



A 



ffiB. 



V 2 



1 1 . e r°°[K^)- 



sin 



eo' 



iTi/q5-i(cj>- f j.TB) cos 2 f -| cot § 



p o 

da—iT J s ^fq^da—i^ 



J° ii$-\£TB) cos 2 |d<7+i^ gin © + 



(38) 



cos 



e 



Jo °° |(^-^i)+^ v ^-i(</»+^TB)sm 2 |-ftanf 



da+iT I \fq2da- 



J" i(<j>+(iTB) cos 2 f *H-i0o ® 



where @o = 0(0), </>o = 0(0) and the infinite integrals in the above formulae are finite. 

Comparing now the formula (|38|) with (|12() and 1)13(1 we see that we have to put D = in 
the formula (|33|) while for the coefficient A we get 



A 



flB-Y ■ ®0 fl^ |(^-^i)-*T^-i(0- M TB)cos 2 §-§cot§ 

n—< 



rO 

da+iT J s ^/q2da-\-i(po 



(39) 



Consequently it is the solution 6 2 I ( s > ^) which will be continued canonically to the sectors 
n + 1 and n + 1 from which it is subsequently continued to +00 of the real s-axis. According 
to the figure this canonical continuation can be done by representing 6 2 i(s,T) first as a 
linear combination of the next two fundamental solutions 62,2(5, T) and &2 2"( s )^) defined in 
the respective sectors 2 and 2. Next the latter two solutions have to be expressed in the 
same way by a pair of fundamental solutions of the sectors 3 and 3 and so on up to the 
moment when the fundamental solutions of the sectors n + 1 and n + 1 enter this procedure. 
Representing the corresponding fundamental solutions in the form 



b 2 , k (s, T) = ft* T)e Tj ^ ^^^(s, T) 
b 2 - k (s, T) = q ~Hs, T)e- lTj ^ ^%- k (s, T) 



(40) 



k = 2,3, 



,n + 1 



this chain of operations can be handled by the following multiplications of matrices 

M = .M..M> ...M n 



Mi 



1 





lX2,l->2 X2,1^2 
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M k 



p0k v 

- ia k,k e ^X 2 - k - 



>k+l 



el3 ~ k X 2 ,k^k+l 



(41) 



iT P fe y/q 2 {s,T)ds 



0k+i=iT [ Jq 2 (s,T)ds, = -iT \ ^ \/q 2 (s,T)ds 



so that 



k = 1, . . . , n 



63,1(8, T) = M 21 6 2jn+ i(s ) T) + M 22 6 2 ^p I (s,T) (42) 



Let us note that the phase integrals defining the coefficients fit and (3%, k = 2, ... n, are 
purely imaginary. Moreover the coefficients aj. k , k = 1, ...,n, become pure real and equal to 
each other while the coefficients (5k become equal to — /3^, k = 2, n in the adiabatic limit 
T — ► +00. To be more precise in these latter statements let s't, k = l,...,n, denote points 
where the antiStokes line emanating from sj,, k = 1, ...,n, crosses the Stokes line passing by 
the points s%, s 2 , s n . Then by the assumption 8°c of Sec. 2 we have 

/ Jq2(s,T)ds = ( 1 Jq 2 {s,T)ds + [ Jq 2 (s,T)ds = 
Js k Js k Js ' k 

j°y q2 ( S ,T)ds+o(±) 

Jq 2 (s,T)ds = r Jq 2 (s,T)ds (43) 

Jsf. Jsi 

I Jq2(s,T)ds = I k+1 Jq 2 (s,T)ds = 

J St. J s' r 



k f s k+l 



k 



+ + K+1 Jq2(s,T)ds = 

'Sk J Sk+l j 

r s k+i / 1 

/ J q2 ( s ,T)ds + 0(-) 

k = 1, n 

i.e. each point sjr, k = 2, n tends to its corresponding limit s^, k = 2, n when T — * +00 
with the rates shown in (|43[). 
Rewriting the Eas.(|33|) as 

a 2 ( S ,T) = Ae U °^- lT " 1 )\ i (s,T) (44) 
and taking into account (|15jl we get 

a_(+oo,T) =M 21 ( — ) sin^x 
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cxp 



f/° 








[5 




[J-oc 







111 B cos cot — 

^ ' 2 2 2 



~ - iTu^j + iT^i - i(<j) + /iTB) sin 2 ® 



«r f° V&da + »T /° y/qjds I (45) 

■/si Js„ J 

since the second term in (|42|) vanishes in the limit s — ► +00 along the real axis (because 

0(+oo,T) = by assumption). 

The formula ([45 1) is just the one which corrects the erroneuos formula (29) of the paper 
(as well as the other formulae corresponding to other cases considered in the cited paper). 
It should be stressed that the formula (|45|) is exact. In this form it looks however very 

complicate because of the complicated structure of the matrix element M^i- The latter is 

polynomial with respect to the coefficients a% k , k = 1, n, and rational with respect to the 

X-coefhcients. Exposing its linear terms in a's we get for it 
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tt _i A *T //* y/^T)ds-iT j°k y/^T)ds+iT J£ y/^T)ds 
J. J. ~X-k—*k / - ^ 



fc=2 



fc=l 



fc-1 



Xfc—*+T II XT— 1+1 II Xz-7+T + 



(46) 



z=i 



l=k+l 



Nevertheless in the adiabatic limit T — > +00 the formulae 1)45(1 and (|46|) simplify greatly 
since then all %'s coefficients of M21 become equal to 1 and in the multiplication of the limiting 
matrices Mjt's, fc = 1, . . . , n, all terms containing powers of the factors k , k = 1, n higher 
than the first ones have to be neglected. Further everywhere where B(s, T) 7^ the adiabatic 
limits T — > +00 of ^ ^~ — iTu?j iiiT-^q^s, T) are exactly the same in their forms as those for 
s — ► ±00 given by (jHZJ). The condition B(s,T) ^ is satisfied obviously for the integrals in 
the formula ()45|) taken along the real axis. However the phase integrals defining the element 
M\2 are taken between zeros of yq^s, T) which in the adiabatic limit coincide with the ones 
of Bq(s). Nevertheless this trouble can be easily avoided by representing the corresponding 
integrations linking pairs of zeros (s^, Sfc), (sj, s^) and (s&, s n ), fe = 1, n, by the ones along 
closed contours C s - Sk , C SlS - and C SfeSn , surrounding respective pairs of zeros. The same idea 
applies to the two integrations between the pair of points (si,0) and (s n , 0) except that the 
corresponding contours C Sl o and C Sn o are not closed but starts and ends at s = points lying 
on two different sheets of R^. Therefore making first use of the Eas. ()43|) we can apply the 
asymptotics p7|) to all the phase integrals in the formulae 1)45 Jl and (jj^j) so that the former 
takes the following form when the integrations along the real axis is performed 



a° d (+oo,T) = -it&n—e 



9o H/c Si!1 + 4// Cs „o W» TB 



e 

sin 



■itj> COS © ds 



fc=i 



l(jc 3l s h Ic SkSn )\ l 



ipTB- 



e 

sine 



■ j<A cos ] ds 



(47) 



where it is assumed that all the quantities are now determined by the asymptotic field Bp. 
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In the present form of the formula (|47|) only the integrations of the middle term — in 
the exponents can be performed explicitly (since / ^^ds = In tan ®). First let us note that 

because Jfe = \ ( glg^ ~ b° +b° 0> 1 ) ' so that lntan § = 2 ln 5+5^ ' we can always choose 
all the integration contours in (|47|) in such a way to avoid possible roots of Bq^ z ± Bq = so 
that the unique singularities which remain inside these contours are branch points of Bq(s). 
Therefore to the corresponding integrals of — along the closed contours C SlSk , C SkSn , k = 
1, ...,n, and C SIS1 can contribute only roots of the function F(s) = ^°|^ + ^°' z |^ . Net results 
of these contributions depends however on detailes of mapping of the s-Riemann surface on 
the F-one. If after such a mapping a closed contour C 7 rounds the zero point of the -F-plane 
n 7 times (we take n 7 to be positive for anticlock orientation of the contour and negative for 
the opposite case) then a contribution of this zero point to the corresponding contour integral 
of — is —mn^. The remaining two open integrals along the contours C Si o and C Sn o can 
contribute only by their limits giving 

= -lntan— + il- (48) 

with some integer I since Bo )Z (s) is regular at the points si,s n and values of Bq(s) on both 
the sheets differ by sign so that F 2 (0) = Ff x (0) where Fi (2 (0) are values (both real) of F(s) 
at s = on the 'first' and 'second' sheets respectively. 
Therefore we obtain the following final result 

JL i . , , -l( f - f ) (iaT B-i6 cos e)ds 

eI »(n SlSfc -n SfcS J e H Jc n s fc J Cs k s n Jy^ v > (49) 

fc=l 

Since cos G = -g^ B o,xB^,y^Bo, y B ,x wg can ^he integrations in (|4^|) to paths 

linking respective points to get 

a a _ d (+oo,T) = _ji+l e -4 ,TO «in e 2 V Js i Ja i v ^ v ; x 

y"e4* 7r ^ a x"fc-' i 'fc»n) e 2 V Js i Js W v v ; (50) 



k=l 



We should remember that all the integrations in (|50[) run along paths avoiding roots of 
the equations Bq >z ± Bq = 0. 

For the corresponding transition probability we obtain 



,ad^ _ -^C(nTB-4>cose)ds 



p« a (T) = e Js i 



n 



y ns l s k ns k s n,l p \ Js l Js k J v ' 



fc=l 



2 



(51) 



Formulae similar to (|50[) and (|51j) have been found by Joye, Mileti and Pfister 5 . In fact 
if we apply the assumptions made in the last paper by its authors these formulae become 
identical, up to an overall phase in Q5U|) . with the corresponding ones found by the authors 
mentioned. 
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The last formulae take on particularly simple forms for the case of two turning points 
lying on the upper Stokes line drawn on Fig. 3. when the equations -Bo, 2 ± Bq = have no 
solutions inside the strip bounded by the two Stokes lines on Fig. 3 and on the lines themselves. 
We can then deform all integration paths in the formula (|50[) to ones along the Stokes and 
antyStokes lines so that the corresponding integrals will have explicitly pure real or pure 
imaginary values. We get for this case 



a a _ d (+oo,T) = -/+ 1 e -l , ™ s i s ie 2 l Js i J n J «a^ " 
,. jonn slS2 -\ HI (v»TB-i* cos 0)ds g - J 



30 ~\ HI (ipTB-ii cos e)ds + e _ 1 innai a2 e | (i„TB-i4 cos &)ds 
.1+1 -jiirn sisi + ±m(+ J°_+ f°)(nTB-<j>coae)ds -|9 J^ 1 (fj,TB-<j)cos e)ds 

cos J (fiTB — <p cos O^j ds — -irn SlS2 ^J (52) 
so that for the corresponding transition probability we get 

s — Q [ sl (uTB-<j>cos0)ds o f I f s ' 2 / ■ \ 1 \ , s 

P ad {+oo,T) = e h i y ' cos 2 ^-Ky (jiTB-^cosQjds--im slS2 j (53) 

7 Examples of NED systems 

An example of a class of fields B with the NED property has been considered recently 
by Berman et al |17] . The fields are defined by putting B z (sT,T) = B^, B x (sT,T) = 
f(s) cos(loqsT), B y (sT,T) = f(s) sin^osT) with /(s) having the properties 1° — 3° of the 
field B and vanishing at both infinities of the real axis. This problem is however unitary 
equivalent to the one with the field B = [/(s),0, B M — ^] so that for this case we have 

B = y / '(f) 2 + / 2 (s) where fl = fiB^ - lo and <p = 0. 

Assuming for /(s) properties desired by the assumption 1° — 8° of Sec. 2 we get using the 
formula (j5U|) 



£ e iMW^ " C) V / (^) 2 ^^ (54) 
fe=l 

where sj and s&, fc = 1, ...,n, are roots of the equations f(s) = 

If there are only two turning points Si and S2 then according to formula (|52[) we get 



ai d (+cx>,T) = -2i m e 4 ' 2 """ V J *i ' J **J V w x 



cos ( ^™£ 2 y(^)" + / 2 (*)^ " \™ aiaa J (55) 



It is now not difficult to establish that the contour C SlS2 rounds the zero point on the 
-F-plane twice (see Fig. 4 and Fig. 5). Therefore we obtain finnally for this case 



a a _ d (+oo,T) = -2i l+ 



1 "i— ,si+¥^(+S°_ + S° 2 )y/(^ 2 + fHs)ds 
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^ J,iVW sm\-fimj M-J +f 2 (s)ds 



(56) 



and for the corresponding transition amplitude 



(57) 



The last two formulae have been obtained earlier by Nikitin and Umanskii ^Sj as well as 
by Crothers and by Davies and Pechukas |5J using the steepest-descent methods. 
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Fig. 4 The integration contour C S1S2 mapped into the /-plane 
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Fig. 5 The integration contour C SlS2 mapped into the F-plane 



As a second example we shall consider again the Nikitin Hamiltonian for the atom- atom 
scattering. The model of Nikitin JH] describes the scattering A*+B^A+B+Ae of the exited 
atom A* moving with a small velocity v with the impact parameter b' and scattered by the 
atom B. The interaction between the atoms is of the dipol-dipol type. The latter example 
was analyzed in the context of the adiabatic limit v — > also by Joye et al [2]. 

The Hamiltonian for this system reads ([IB|> paragraph 9.3.2 and |19j): 



H(R) 



' Ae 


c 


2 


W 


C 


Ae 


. w 


2 



(58) 
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\Jy 2 + v 2 t 2 is the distance between the atoms. 
Introducing d = (2C/Ae)t as a natural distant unit for this case and T = d/v as the 
corresponding adiabatic parameter and rescaling: t — ► sT and 6' — > bd we get from (|58|1: 



where Ae and C are constants and R 

i 
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H{s) = — 



Ae 
2~ 



1 



(£> 2 +s 2 )2 



(fc 2 +s 2 )2 



(59) 



(& 2 + s 2 ) 



0,1 



so 



(0,0, l)f . 



In the 'magnetic field' language we have of course B(sT, T) - 

that all the assumptions 1° — 8° above are satisfied with B =t: (T) = B =t (±oo, T) 

Obviously the last form of the B-field shows that it belongs to the class of Berman et al 
with two turning points on the "main" Stokes line (see Fig.l) so that the formulae ()56|) and 
(foTj) are applicable readily. 



8 Discussion and conclusions 

In our present calculations of the adiabatic limit for the transition amplitudes in the 
two energy level systems we have corrected erroneous formulae of our previous paper pQ. We 
have considered systems with the NED properties, i.e. for which their corresponding Stokes 
graphs do not differ essentially from their adiabatic limit forms. 

A formula (|5U|) which gives the corresponding transition amplitudes in the adiabatic limit 
shows that these amplitudes result as an interference of contributions coming from all complex 
conjugated pairs of turning points lying on the same complex conjugated Stokes lines of the 
respective limiting Stokes graph. Up to an overall phase it coincides with the one of Joye, 
Mileti and Pfister 0. 

A particularly simple formula for the transition amplitudes follows from a general one 
(|52j) when the latter is applied to the NED systems considered by Berman et al |17j with 
two turning points on the "main" Stokes line. Namely, it obtains then the form ()56j) found 
earlier by Nikitin and Umanskii JH] as well as by Crothers jSHj and by Davies and Pechukas 
[2*T] using the steepest-descent methods. 
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